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ABSTRACT 

In gravitostatics, the minimally generalized Newtonian theory, the simplest cosmological model implies a 
cosmological redshift with acceleration parameter —1. If the gravitational potential satisfies a wave equa- 
tion, the total mass density (visible plus dark) is one third of the critical density, in good agreement with 
observation. There is no horizon problem and no need to invoke dark energy. 

Subject headings: gravitation — cosmology: theory — galaxies: high redshift — dark matter — early universe 



1. Introduction 

The empirical evidence that might distinguish between different, relativistic theories of gravity is very sparse. 
Any acceptable theory must of course be compatible with the Newtonian theory. It used to be thought that 
ON \ the small deviations from the Newtonian predictions, the so-called 'classical tests of relativistic gravity', can 

be used to discriminate between relativistic theories. However, a slightly modified Newtonian theory, called 
gravitostatics, accounts for these deviations. Any relativistic theory must satisfy the classical tests, but this 
shows only that it is not incompatible with gravitostatics. (One must emphasize that gravitostatics is not a 
relativistic theory, although it reduces to special relativity in the limit when all gravitational fields vanish.) 
CO . To choose between different compatible theories, other experiments are needed. 

CO 

In addition to the classical tests, there are other, small, 'postnewtonian'deviations from the Newtonian 
predictions, but they yield only weak constraints on possible theories. The most promising testing ground 
for relativistic theories is the study of gravitational radiation. Although experiments are very difficult, there 
is one well-established, observational result: the rate of energy loss by gravitational radiation from a binary 
pulsar. This agrees with the prediction of the Einstein theory, and is the only good piece of evidence for it 
at present. The Einstein theory is widely accepted for aesthetic reasons, and because of the misapprehension 
»! . that the classical tests support it uniquely, as discussed in the last paragraph. 

Cosmology is another field for the testing of gravitational theories. The Einstein theory had difficulties here. 
There was the 'horizon problem': spatially separated regions are causally disconnected at early times. There 
was also the 'flatness problem': unless the initial conditions are chosen with extreme precision, the universe 
does not unfold as it should. Both difficulties were ingeniously overcome by the theory of inflation. It is still 
unclear whether this produces exactly the correct conditions for the subsequent expansion of the universe; 
even if it does, many feel that it is too contrived and that it destroys the aesthetic appeal of the theory. 
Recently, measurements of very distant objects have shown that the Einstein theory is admissible only if the 
universe is filled with an otherwise unobservable 'dark energy'. The invention of such entities is a mark of 
theories under strain: one thinks of phlogiston and the ether. 

In this uncertain situation, it is illuminating to see what gravitostatics says about cosmology. This requires 
another small generalization of the theory, the Poisson equation, inherited from the Newtonian theory, must 
be modified by the addition of a term involving time derivatives. In this paper, the cosmological predictions 
of the modified theory are shown to be free from the problems of the conventional theory. The way is now 
open for the investigation of properly relativistic theories with the same qualities. 

The rationale for this paper is that, since the empirical evidence for the Einstein theory is not strong, it 
is worthwhile to investigate other theories. This will reveal a wider range of possibilities and will help to 
distinguish real features of the world from merely theoretical constructs. It is suggested that inflation and 
dark energy may be examples of such constructs. The opposing view, more widely held, is that there is 
no point in developing other theories so long as observations do not clearly contradict the Einstein theory. 
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Which opinion one adopts is a matter of taste or temperament. It requires cool assessment, not religious 
fervor. 

2. Gravitostatics 

It should be well known that the Newtonian theory of gravitation can be generalized to be compatible with 
special relativity in the limit of vanishing gravitational fields, and to account for the gravitational redshift 
and the other classical tests of relativistic gravity. The generalized theory, gravitostatics, is non-speculative, 
in that it makes no new hypotheses but simply takes the old ones more seriously. For example, one finds 
that the assumption that the gravitational potential is arbitrary to the extent of an additive constant implies 
that the rate of a clock depends exponentially on the gravitational potential — and that the units of length 
and mass have a similar dependence. A full account of all this can be found in (Rastall 1991); a shorter, 
more accessible source may be (Rastall 1975). The next paragraphs try to summarize enough of the theory 
to make the rest of the paper comprehensible. 

In gravitostatics, there are preferred coordinate systems, called Newtonian frames, which are similar to the 
Galilean frames of Newtonian theory or the inertial charts of special relativity. As in Newtonian theory, the 
gravitational field is described by a function $, the gravitational potential. The measured values of universal 
constants such as the speed of light c the Newtonian gravitational constant G, and Planck's constant h are 
assumed to be independent of $. In a static gravitational potential (that is, when $ may depend on the 
space coordinates x = (x , x 2 , x 3 ) of a Newtonian frame but not on the time coordinate t) the time interval 
measured by a fixed standard clock is not in general equal to the change in the time coordinate t. Similarly, 
the distance as measured with a standard measuring rod between the spatial points with coordinates x and 
y is not in general |x — y| = [(x m — y m )(x m — y" 1 )} 1 ^ 2 (the summation convention is assumed unless stated 
otherwise: sum the repeated index to over its range (1, 2, 3)). The spatial geometry as measured by standard 
measuring rods is consequently not Euclidean. Directly measured times and lengths will be termed empirical; 
those defined in terms of the space and time coordinates of a Newtonian chart will be termed Newtonian. 

The energy of an object at any point can be measured in terms of the energy of a reference object at a nearby 
point (the reference object might be a standard slug of metal or a specified atom in its ground state, for 
example). This is called the empirical energy. From the laws of mechanics, which are a trivial generalization 
of the special relativistic laws, it follows that there is an energy function which is constant along the path of 
a particle in a static gravitational potential. This is called the Newtonian energy, and is in general different 
from the empirical energy. 

Empirical quantities may be denoted by a subscript E and Newtonian quantities by a subscript N. For local 
quantities — those which are defined in a region where $ is almost constant — one proves that the ratio of 
a Newtonian quantity to the corresponding empirical quantity is an exponential function of $ (this follows 
from the assumption that $ is arbitrary to the extent of an additive constant). If tjy, £n, &n are a time 
interval, length, and energy measured in Newtonian units, and tg, £e, &e are the corresponding quantities 
in empirical units, one has 

t N = t E exp[-r((j) - (p )], £n = £e cxp[— A(</> — </ > o)] > (1-1) 
e N = e E cxp[-rj((f) - </> )], 

where r, A, n, and 0o are constants, 4> = <&/ce 2 , and ce is the speed of light in empirical units. One calls <j> 
the dimensionless potential. Newtonian and empirical values are equal when <p — <Po\ a Newtonian frame for 
which this is true is called a (f>o frame. 

The magnitude of a dimensionless quantity is independent of the units of measurement. It follows that, 
when measured in Newtonian units, the magnitudes of two quantities with the same dimensions must have 
the same dependence on (f> (because their ratio is dimensionless). Since the dimensions of a mass m are 
[to] = [e£~H 2 ], eq.(l.l) gives 

m N = m E exp[-ij,((p - fo)], (1.2) 
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where /i = rj — 2A + It. To avoid an excess of subscripts, one makes the convention that universal constants 
without a subscript E or N are to be interpreted as empirical quantities (the speed of light c is to be 
interpreted as Ce, Planck's constant h as He, etc.). Anything else written without a subscript E or N (i.e. 
anything but a universal constant) is to be interpreted as a Newtonian quantity, unless stated otherwise (the 
velocity V of a particle is to be interpreted as Vjv, etc.). 

We have now set up the conceptual framework of gravitostatics. To go further, we must learn how to 
calculate the motion of a particle in a gravitational field. This would seem to require some additional, 
arbitrary hypothesis. However, the Lagrangian method allows us to proceed in a non-arbitrary, or at least 
very natural way. 

The Lagrangian of a free particle in special relativity is L = —mc 2 /-f, where m is the constant and invariant 
proper mass, 7 = (1 — X /2 c -2 ) -1 / 2 , and V = |V| is the speed. The guiding principle of gravitostatics is to 
adopt the laws of special relativity and Newtonian gravity with minimal changes. One therefore assumes 
that the Lagrangian of a particle subject only to a gravitational field is the same as that for a free particle 
— except that everything is to be measured in Newtonian units. With the previous conventions, one has 

L = -mc N 2 h, 1 ={l-V 2 c N - 2 )- 1 ' 2 , (1.3) 

where m is the proper mass in Newtonian units. As in Newtonian mechanics, it is assumed that tue, the 
proper mass in empirical units, is constant along the path. (It would of course be more correct to write the 
Lagrange function as L(Z, V, •), where Z is the position of the particle, and its value at time t as L(Z(t), 
V(t),t).) From (1.1), one gets 

cjv = cexp[(r — A)t/>], mcjv 2 = thec 2 exp(—r/ip), (1-4) 

where ip — <f> — 4>o- If ip and V are small, the Lagrangian must be equivalent to the Lagrangian m,EV 2 /2 — 
rriEC 2 <j> of Newtonian gravitation, which implies that 77 = —1. 

The r component of the momentum of the particle is p r = dL/dV r = m-fV r and the Euler-Lagrangc equation 
is dpr/dt = dL/dZ r . Since V r p r — L = m-f(V 2 — cn 2 J~ 2 ) = mjCN 2 , one has 

{d/dt){m 1CN 2 ) = -8 t L, (1.5) 

where d t L denotes the partial derivative with respect to the last variable in L = L(Z, V, •). If <f) is static, 
then d t L = 0, and E = tujcn 2 is constant along the path of the particle. It is identified with the energy 
of the particle in Newtonian units. The Newtonian unit of energy can thereby be defined everywhere in the 
Newtonian frame. One may choose the Newtonian unit to be the same as the empirical unit in a region 
where the potential is fa. (If 4> is quasistatic, as will later be allowed, then (1.5) gives the change in E along 
the path.) 

Rapidly moving objects are permissible in gravitostatics, provided that they do not contribute significantly 
to the gravitational field. One can introduce a classical model of a photon as the limiting case of a particle 
whose mass tends to zero and whose speed tends to the speed of light in such a manner that its energy 
E = mjCN 2 remains bounded. The momentum of the photon is parallel to its velocity and has magnitude 
E/cn- In a static potential, the Newtonian energy E of the photon is constant along its path. 

A photon of Newtonian energy E is associated with an electromagnetic wave of Newtonian frequency v. 
Planck's law states that Ee/ve — h, where h is Planck's constant, or equivalently that E/v = h^. In 
a static potential, the Newtonian energy E of the photon is constant along its path, and the Newtonian 
frequency v of the associated wave is constant (otherwise, the number of wavelengths between two spatial 
points is not constant in time). Hence E/v is constant along the path. At the potential <j>o, where Newtonian 
units are the same as empirical, one has E/v = h and it follows that, at any potential, E/v = = h: the 
value of Planck's constant is the same in Newtonian and empirical units. Since /ijv = /iexp[— (ry + r)^], from 
(1.1), we have r = — rj = 1. 
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In the Newtonian theory, the gravitational potential $ obeys Poisson's equation d m d m <& = 4irGp, where G 
is the gravitational constant and p is the mass density of the sources. In gravitostatics, the only difference 
is that one uses the energy density e of the sources: 

d m dm4> = 4:irk N e, (1.6) 

where (f> — <I>c~ 2 and fcjy = Gncn~ 4 , or k = Gc~ 4 (note that $ is measured in empirical units). 

In Newtonian theory, the active gravitational mass of an object, which determines the strength of the 
gravitational field which it produces at large distances, is the same as the passive gravitational mass (the 
weight) and the inertial mass. In gravitostatics, similarly, the active gravitational energy is assumed to 
be equal to the total energy of the sources. A short calculation shows that this implies A = 77. All three 
constants in (1.1) arc now determined: 

r = -A = -77 = l, (1.7) 

and eqs.(1.4) and (1.3) become 

c N = ce 2 ^, mc N 2 = m B cV, L = -m B cV(l - V 2 cr 2 e^f 12 . (1.8) 

Using the values (1.7), one calculates the gravitational redshift, the bending of light by the Sun (and the 
related radar time-delay), and the anomalous perihelion shift of Mercury. These are often called the classical 
tests of relativistic gravity. In fact, they are predicted by gravitostatics, which is only a slight generalization 
of the Newtonian theory. 



3. The cosmological redshift 

Like Newtonian theory, gravitostatics can be applied when the gravitational field is slowly varying in time. 
The simplest cosmological models are those in which the gravitational potential is independent of the space 
coordinates and is a function only of the time coordinate t of a Newtonian frame. We take the special case 
when the potential is a linear function of the empirical time T: 

iP = <f>-<j> = -HT, (2.1) 

where (f>o and H > are constants and T is a function of t only. Since T at any spatial point can be regarded 
as time measured by a standard clock fixed at that point, one has dT/dt = = e~ HT from (1.1) with r = 1. 
Integrating and choosing T = when t = 0, one gets e - ^ = Ht + 1. As t — ► 00, tp — ► —00 and T — > 00; 
as t — ► —\/H, tp — ► 00 and T — ► —00. In this model the 'horizon problem' of cosmology does not arise: 
since cjv = ce 2 ^ , signals can, in principle, be sent to any spatial point x at a time t from any other spatial 
point x' at a time t', provided that t' is sufficiently close to the value —1/H. 

To define the redshift, we suppose that light emitted by an atom at rest at the spatial point 1 in a Newtonian 
chart has frequency ve, measured in empirical units. It travels to the spatial point 2, where its frequency is 
measured to be vie- This is compared with the frequency ve of the light emitted from an identical atom at 
rest at the spatial point 2. The redshift z (of the light from 1 as measured at 2) is defined by 

z = [y E - v 2E )/v 2 e = (E E - E 2 e)/E 2 e, (2.2) 

where the frequency ve of the light and the energy Ee of the corresponding photon are related by ve = Ee/H, 
and similarly v 2 e = E 2 e/Ii. 

It is easiest to calculate the redshift by again considering the motion of a classical 'photon'. The Newtonian 
energy of such a photon changes according to (1.5). From (1.8) one has dL/dip = —E(l + V 2 cjy~ 2 ), and 
putting V = cn gives 

dE/dt = 2Ed t ip. (2.3) 



4 



As before, a photon of empirical energy Ee is emitted by an atom at rest at the spatial point 1 at the instant 
ti in a Newtonian chart. The Newtonian energy of the photon is E\ — Ee^ 1 , from eqs.(l.l) and (1.7), 
where ipi = (f>i — (fio an d (po is a constant. The photon travels to the point 2, where its Newtonian energy is 
E 2 . From (3.3), \n{E 2 /Ei) = J 2d t ip(Z(t),t)dt, where the integral is from tito t 2 . The photon's empirical 
energy at point 2 is 

E 2E = Ehe-** = £ B exp[^i - ip 2 + [' 2d t ip(Z(t),t)dt}. (2.4) 
From (3.2), the rcdshift is 

z = (E E /E 2E ) - 1 = cxp[^ 2 - Vi - / 2 2$V(Z(f)> *)<**] - 1- (2-5) 

In the special case of a static potential, one has d t ip — and the redshift is z — exp(?/>2 — tpi) — 1 w ^2 — V'l 
if |^2 — ^1 1 -C 1. Note that the redshift is positive if the light goes from a lower to a higher potential. If 
d m ip = everywhere, so that ip is a function of i only, eq.(2.5) becomes 

z = cxp(V>i - ip 2 ) - 1- (2.6) 

Here the redshift is positive if the light goes from a higher to a lower potential. 

In the special case when ip = —HT, eq(2.1), the redshift is 

z = cxp[H(T 2 - T x )] - 1 = H(T 2 - Ti) + (l/2)H 2 (T 2 - Ti) 2 . (2.7) 

The last equation is valid if H(T 2 — Ti) is small. The term in H 2 corresponds to a deceleration parameter 
q = — 1, in the conventional terminology (see eq. (6.2) or eq.(29.15) of (Misner, Thornc, & Wheeler 1973)). 
Eq.(2.7) implies that the universe is 'accelerating' — it expands more rapidly at later times. This agrees 
with recent measurements on the redshifts of distant supernovae and galaxies (Knop, et al. 2003), (Tonry, 
et al. 2003), (Daly & Djorgovski 2003). 

4. A field equation and dark mass 

The field equation (1.6) is the same as the classical Poisson equation, except that the mass density of 
the sources is replaced by their energy density. An apparent defect of the equation is that the energy 
density of the gravitational field does not appear in the source term. It is however easy to remedy this: 
instead of choosing the field variable to be 0, or ip — <j> — cf>o, one takes it to be e^/ 2 . One then has 
2e^/ 2 d m d rn e^/ 2 = d m d m ip + {l/2)d m ipd m ip and (1.6) becomes 

2e-^l 2 d m d m e^' 2 = 4irk N (e + e G ), (4.1) 

where ec — (87rfc) dmipdmip is the energy density of the gravitational field. Since the dimensions of 
k = Gc are length divided by energy, we have = k, from (1.1) with rj = X = — 1. Note that ea is never 
negative — unlike the energy density in Newtonian gravitation, which is —eg. 

Eq. (4.1) is completely equivalent to (1.6). It is the equation satisfied by a static potential ip. If the potential 
is not static, then the field equation must involve time derivatives if it is to determine the time evolution of 
the potential. Here a speculative clement enters the theory. We guess, as the most obvious possibility, an 
equation similar to the wave equation: 

2e^/ 2 [d m d m e^' 2 - cn-^Cn- 1 ^ 2 )] - 4^fc(e + e G ). (4.2) 

(The cjv are inserted to ensure that the equation is invariant under change of <po chart.) 



5 



If we assume that ip(t) — —HT, as in (3.1), and use cm = ce 



^ and dT/dt = e*, we have 



2e-^/ 2 c N - 1 d t (c N - 1 d t e^/ 2 ) = -(H 2 /2c 2 )e^. 



Defining the critical density e c by e C E = 3H 2 /8irkc 2 , and e c = e C Ee we find that (4.2) becomes 



Since the space derivatives vanish in this cosmological model, we have ec = and e = e c /3. 

The average density of the visible matter in the universe is only a few percent of the critical mass density. 
Observations of the motions of galaxies and of gravitational lensing effects, etc., show that 'dark' gravitating 
matter must be present. The average total density of visible plus dark matter is about 30 percent of the 
critical mass density. If we identify e = e c /3 as the average energy density of all the matter, visible plus 
dark, we have quite good agreement with observation. The theory does not say anything about the nature 
of the dark mass. 

5. Conclusions 

As emphasized earlier, gravitostatics is not a new, speculative theory, but only a minimally modified form 
of Newtonian gravitation. In this paper, we have generalized it slightly by allowing gravitational fields to 
vary slowly in time, and have assumed that they satisfy a wave equation. The simplest solution of this 
equation implies that the average total mass density is one-third of the critical mass density, in agreement 
with observation. The universe has an infinite past, as measured in empirical units, and there is no horizon 
problem. The deceleration parameter has the value — 1, which implies that the expansion rate of the universe 
increases with time, and which accords with observations of objects with high redshift. One does not have 
to postulate a mysterious 'dark energy'. 

An obvious generalization of the results is given in the following Appendix. We showed above that if 
ip = —HT then the time-dependent term in the field equation (4.2) gives rise to a contribution to the energy 
density of the sources which is constant in time when measured in empirical units. If we now write the time- 
dependent term as 2e~^/ 2 c N ~ 1 dt{cN~ 1 d t e^/ 2 ) = cr 2 e~ 2 ^B, our result is that B is constant if tp — -HT. 
One naturally asks what is the most general tp that implies a constant B. We find that, in the general case, 
the contribution to the energy density may be less than before. The deceleration parameter may depend on 
the time of observation — it lies between —1/2 and —1 and tends to the latter value at large times. 

This paper uses the conventional model of a universe with homogeneous space sections, which cannot be 
expected to accord exactly with obervation. The equations of any such model will differ from those found by 
averaging over our actual, lumpy universe, because the non-linear terms will almost certainly not average to 
zero. Once the distribution of dark matter is known, it should be possible to estimate the size of the error. 

In summary, we have shown that a simple, rather primitive theory of gravitation can account for the large- 
scale structure of the universe without invoking hypotheses such as inflation or dark energy. The challenge 
now, for vigorous young theorists, is to find a properly relativistic theory that is at least as successful. This 
is not trivial because one cannot assume that the potential is a relativistic invariant (it is well known that 
relativistic scalar theories of gravity do not work). The examples of exact theories in Chapter IX of (Rastall 
1991) do not take account of the field equation (4.2). 



We shall find the most general function %p that depends only on the time coordinate t and that satisfies 
the equation 2e~^/ 2 dt(e~ 3 ^/ 2 d t ?p) — B, where B is a constant. We write w — e -3 ^/ 2 and the equation 
becomes d 2 w/dt 2 — — (3B/2)w~ 1 / 3 . Multiplying by dw/dt and integrating gives dw/dt = ±(a + bw 2 / 3 ) 1 / 2 , 



2e-^' 2 d m d m e^l 2 = 4irk(e + e G - e c /3). 



(4.3) 



6. Appendix 
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where b = —9B/2 and a is a constant of integration. The upper sign is chosen so that dip/dt < 0. Defining 
u = w~ 2 / 3 = e^, gives t = —(3/2) J u~ 2 (au + b)^ x / 2 du, which is a standard integral. 

The proper time T satisfies dT/dt = e^. Since dt/du = -3/2u 2 5, where S = (ae^ + b) 1 ' 2 , and du/dip = u, 
we have dT/dip = -3/25, and hence V' = d^jdT = -25/3, ip" = (2/9)ae^. Eq.(2.6) implies that the 
redshift observed at time T 2 is 

z = -ip'(T 2 - Ti) + (l/2)(^' 2 + V")(r 2 - T,) 2 , (6.1) 

provided that T 2 — T\ is small, where ip' and are to be evaluated at T 2 . The conventional form of this 
redshift is 

z = H(T 2 - + (1 + q/2)H 2 (T 2 - T^ 2 , (6.2) 

where H and q are the Hubble parameter and deceleration parameter at T 2 . It follows that H = 25/3 and 
q = —1/2 — b/2S 2 . Note that q lies between —1/2 and —1, and that q tends to —1 when T 2 is large. If 
T 2 = T , the present epoch, and if the Newtonian frame is chosen so that tp(T ) = 0, then Newtonian and 
empirical units coincide at T , and S(T ) = (a + b) 1 / 2 . Writing H = H and q = q Q , we have 

a=(9/2)# (l + «b), b=-(9/4)(l + 2q )H . (6.3) 

This determines a and b in terms of the measurable Ho and go- the earlier assumption that ip = —HT, 
eq.(2.1), corresponds to go = — 1- 
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